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^ ■ Abstract 

Detection of 5 events by the Liquid Scintillation Detector (LSD) on February, 
' 23, 1987 was interpreted in the literature as the detection of neutrinos from the 

first stage of the two-stage supernova collapse. We pose rigid constraints on the 
properties of the first stage of the collapse, taking into account neutrino flavour 
' conversion due to the MSW-effect and general properties of supernova neutrino 

emission. The constraints depend on the unknown neutrino mass hierarchy and 
mixing angle ^13. 
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O 1 Introduction 



SN1987A was the only supernova to date which produced a measured neutrino signal. 



' Four experiments reported the detection of neutrinos: LSD (Liquid Scintillator Detec- 

^ ' tor, [1],[2]), KII (Kamiokande II, [3]), 1MB (Irvine-Michigan-Brookheaven,[4J) and EST 

(Baksan Scintillator Teleskope, p]). While LSD registered neutrino burst at 2:52UT, 
February 23, the other three experiments - at 7:35UT February 23 (UT stands for Uni- 
tary Time). Each experiment reported only one burst: LSD observed no statistically 
significant counterparts for the KII, 1MB and BST neutrino signals and vice versa. This 
puzzhng discrepancy could be, in principle, explained by a two-stage supernova collapse 
hypothesis, as was stated in [Sj-|Tn]. Various two-stage supernova collapse models, pro- 
posed in this papers, implied the transition from a protoneutron star to a black hole 
and/or the formation and evolution of a close binary system inside the exploding star. 
Properties of the LSD signal and non-observation of its counterpart in other detectors 
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Table 1: Type, material and mass (in tons) of the detectors, numbers of events (Ng^) 
at 2:52 and 7:35 (according to the cited references). Only events with measured energy 
exceeding the energy thresholds of corresponding detectors were taken into account. 





LSD 


KII 


1MB 


Baksan 


Type 


scintillator 


cherenkov 


cherenkov 


scintillator 


Material 


C„H2„ 90t 


H2O 2140t 


H2O 5000t 


C„H2„ 200t 




Fe 


200t 








iVev at 2:52 


5 


1, 2J 


2 [SIE] 


m 


1[5J 


A^ev at 7:35 


2 yioj 


11 m 


8 


615] 



allow to pose rigid constraints on the properties of the first stage of two-stage collapse 
scenarios, MSW effect [HI [12] in the matter of the star being of crucial importance. For 
example, it was shown [13] that accounting for the neutrino flavour conversion can spoil 
the reported [lOj concordance of the rotating collapsar model with the data. The analysis 
independent from the particular collapse model was made in [H] . In the present paper we 
extend the analysis of [H]: we conduct a more elaborate statistical study of the data and 
account for supernova shock wave effect, which may influence neutrino flavour conversion 
[TH] . The obtained constraints on the first stage of the two-stage collapse depend on the 
unknown neutrino mass hierarchy and mixing angle ^13. 

The structure of the paper is the following. In Section 2 we describe neutrino signals 
in the detectors involved. In Section 3 general properties of the supernova neutrino fluxes 
and flavour transformations are briefly reviewed. In Section 4 the analysis procedure and 
results are presented. In Section 5 the conclusions are summarised. 

2 SN1987A neutrino signal detection 

2.1 General detector characteristics and SN1987A neutrino sig- 
nals 

Detector characteristics and numbers of registered events at 2:52 UT and 7:35 are given 
in Table 1. LSD reported their event cluster at 2:52 to be a burst, while KII and BST 
regarded their few events at 2:52 as "background fluctuations". On the other hand, KII, 
1MB and BST reported only event clusters at 7:35 to be "bursts". The duration of all 
bursts was about 10 seconds. 

In this paper we discuss only the first stage of the presumable two-stage collapse 
and, accordingly, only the first neutrino signal, which occured at 2:52 UT. Moreover, 
we compare only LSD and KII signals, and do not use 1MB and BST data. Even such 
restricted analysis leads to very rigid constraints on the first stage of the collapse, and 
additional data do not affect the conclusions significantly. 

Time sequence and measured energies of LSD and KII events are listed in Table 2. Also 
the energy threshold Eth and the mean number of background events with energy greater 
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Table 2: Energy thresholds Efh and mean numbers of background events in a ten-second 
interval Ubg according to the cited references; time sequence [t is time in seconds) and 
energies E of events at 2:52 for LSD P[2] and KIT [7J. 





LSD 


KII 


Eth (MeV) 


5 ilJi2j 


7.5 [16] 




0.1 [1J[2J 


0.22 [3][16J 


events 


t (s) E (MeV) 

36.8 7 
40.7 8 
41.0 11 

42.7 7 

43.8 9 


t (s) E (MeV) 

22 12 
30 7.5 



than Eth in a ten-second interval are given for each detector. It should be reminded that 
due to the problems with clocks in KII the timing of KII event sequence could be shifted 
within one minute According to [7] a pair of KII events given in Table 2 is the largest 
cluster of events in KII at 2:52±lmin. 

The imitation rate from the background for the LSD event cluster was fairly small - 
0.7 per year [1],[2]. This justifies the attempts to find an explanation for the LSD neutrino 
signal. 

2.2 Reactions responsible for the neutrino registration in LSD 
and KII 

Neutrinos and antineutrinos can be registered through their interactions with nuclei and 
electrons. Reactions which were responsible for neutrino and antineutrino detection in 
LSD and KII are listed in Table 3. In this Table / = e,/u,r; x = /i,r, and superscript 
"*" denotes the exited states of the nuclei (which immediately decay to ground ones 
emitting nucleons and gammas). We have listed only those reactions which proved to be 
essential in the analysis of KII and LSD signals. Thus we omitted interactions of electron 
antineutrinos with nuclei and electrons because they could at best produce ten times less 
events than reaction (1) (see references in Table 5). Also we omitted all neutral current 
reactions with oxygen because KII, being a Cherenkov detector, was not sensitive to them. 

It should be stressed that while in case of reaction (1) energy measured in the detector 
(i.e. energy of the outcoming positron) almost equals the incident neutrino energy, it is 
not the case for other reactions, especially for those which are accompanied by the nuclear 
excitations. In fact, for LSD it could be exclusively nuclear excitation energy which was 
measured. This issue is discussed in detail in [T0].p!7]. 
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Table 3: Reactions responsible for the neutrino registration in LSD and KII. In the last 
two columns the detector efficiencies are listed and corresponding references are cited. 
They correspond to neutrino energies greater than 30 MeV for reactions (4), (7) and (10), 
and to energies greater than 15 MeV for other reactions. See text for further explanations 
and references. 





reaction 


?](LSD) 


r/(KII) 


(1) 


z/e + p e"*" + n 


1 [IJ 


0.9 ^ 


(2) 


Ue + Fe e + Co* 


0.53 [lOj 


— 


(3) 


i/e + C ^ e + N* 






(4) 


z/e + ^ e + F* 




0.9 [16J 


(5) 


ui + Fe ^ 1^1 + Fe* 


0.53 [lOj 




(6) 


ui + C ^ ui + C* 


1 W 




(7) 


ui + e vi + e 


1 fT] 


>0.5 [16] 


(8) 


Vx + Fe — s> + Fe* 


0.53 [lOj 




(9) 


+ C ^ Px + C* 






(10) 


Ux + e Vx + e 




>0.5 [16J 



Table 4: Electron antineutrino detection effective area Sy^ (in 10 cm^) for KII and 
LSD. 



Ep^ - Arunp - rUe (MeV) 


6 


7 


8 


9 


10 


11 


12 


^p.(LSD) 


0.30 


0.39 


0.50 


0.62 


0.75 


0.89 


1.05 


5p.(KII) 


0.6 


1.8 


3.6 


6.0 


9.1 


11.6 


14.4 



Table 5: Cross-sections for reactions (l)-(6) (in 10 ^*'cm^) according to the cited refer- 
ences. 



£;^(MeV) 


(1) m 


(2) m 


(3) [22] 


(4) m 


(5) m 


(6) p] 


15 


0.162 


0.064 





0.001 


0.021 





20 


0.290 


0.293 


0.003 


0.002 


0.069 


0.003 


25 


0.447 


0.733 


0.019 


0.008 


0.151 


0.012 


30 


0.630 


1.40 


0.050 


0.02 


0.285 


0.027 


35 


0.835 


2.36 


0.095 


0.05 


0.489 


0.048 


40 


1.058 


3.71 


0.151 


0.12 


0.786 


0.073 


45 


1.298 


5.55 


0.218 


0.24 


1.19 


0.101 


50 


1.551 


7.98 


0.292 


0.42 


1.72 


0.131 


55 


1.816 


11.0 


0.370 


0.67 


2.39 


0.162 


60 


2.091 


14.8 


0.452 


1.0 


3.20 


0.195 


65 


2.374 


19.2 


0.530 


1.46 


4.15 


0.224 


70 


2.664 


24.2 


0.608 


2.20 


5.25 


0.254 
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Let us define an effective area S,,^(^p^) for neutrino i/i (or antineutrino ui) detection: 

Suiiui) = T.'qiNiai. 

Here again I = e, fi, r, while i numerates those reactions from the above list which are 
relevant for ui (or detection, r]i is the detector efficiency specific for the z-th reaction, 
Ni is the number of corresponding target particles (electrons, protons or nuclei) in the 
detector and cTj is the corresponding cross section. Effective area depends on energy via 
cross sections and efficiencies. Evidently, mean number of registered neutrinos of a given 
type is proportional to the effective area of the detector. Thus effective area is a measure 
of the detector sensitivity. 

First let us briefly discuss the possibility that LSD signal was due to the flux of electron 
antineutrinos, with energies in the range 6 MeV < Ep^ — ArUnp — rrif. < 12 MeV, which 
corresponds to measured event energies. Here Am^p = 1.3 MeV is the the neutron-proton 
mass difference, and nie = 0.5 MeV - the electron mass. Effective areas for detection 
in LSD and KII through the reaction (1) are given in Table 4. The efficiency for KII, 
which strongly depended on energy in this energy range, may be found in UH]- Efficiency 
for LSD is considered to be 99% according to [2]. The cross section for the reaction 
(1) is taken (here and in what follows) from [18] (eq. (25)). To interpret the data one 
should remember, that energy measured in LSD contained 1 MeV from e~^e~ annihilation 
additionally to positron kinetic energy Ep^ — ArUnp — me- 

Numerous studies (see, for example, [3], [9], [T9], [20]) indicated that if LSD signal 
were due to (7-12) MeV electron antineutrino flux, then 

(1) expected average number of events in KII would be from 5[20j to 28 13], depending on 
the analysis, and 

(2) an enormous energy > lO^^erg, which is greater than binding gravitational energy of 
a neutron star, would be released during the collapse. 

Physical reasons, underlying difficulties in such an interpretation of the LSD data, are 
discussed in detail in [5|^[ I what follows we do not consider this interpretation. 

In 1987 only reaction (1) was regarded to allow supernova neutrino detection. Another 
possibility, on which we will focus our attention, was elaborated only in 2004 [TO] [17]: 
the idea was that neutrinos (not antineutrinos) of sufficiently high (30-50 MeV) energy 
produced the signal in LSD through reactions on iron and carbon nuclei. 

Efficiencies for major reactions relevant for the detection of neutrinos and antineutrinos 
of such energies are hsted in the last two columns of Table 30. Cross sections for reactions 
(1) - (6) in the energy range (15-70) MeV are listed in Table 5. Effective areas of LSD and 
KII for neutrino detection are given in columns 2-5 of Table 6. In columns 6 and 7 one 
can see effective areas for the detection of high energy (with kinetic energy > 14 MeV) 
electron from the z^e^ and elastic scattering in KII. No such electrons were registered 
in KII which should be taken into account in the statistical analysis. The last column 
in Table 6 provides detection area for KII, which is usefuU to constrain z/g flux in the 

^In [5] a two-stage collapse scenario was proposed, which allowed to eliminate this difficulties, but 
implied an extremely heavy (> IOMq) collapsing core. 

^To estimate the efficiency of vi and i>x detection in KII due to the elastic scattering on electrons 
(reactions (7) and (10)) one should take into account the spectrum of scattered electrons. 
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Table 6: Effective areas (in 10^^'' cm^) for neutrino (columns 2-5) and electron antineu- 
trino (column 8) detection, S^^ and Sp^, for KII and LSD; effective areas for the detection 
of high energy (with kinetic energy > 14 MeV) electron from the and elastic 
scattering in KII (columns 6 and 7). See text for further explanations. 



(MeV) 


Ue, LSD 


Ue, KII 


l/y, LSD 


Uy, KII 


HE Uee, KII 


HE z/^e, KII 


l^e,KII 


15 


0.095 


0.99 


0.024 


0.16 


0.060 


0.0085 


28 


20 


0.42 


1.4 


0.087 


0.21 


0.36 


0.052 


51 


25 


1.1 


2.1 


0.21 


0.26 


0.67 


0.097 


79 


30 


2.1 


3.2 


0.42 


0.31 


0.97 


0.14 


114 


35 


3.7 


5.6 


0.72 


0.36 


1.3 


0.19 


155 


40 


5.8 


10.4 


1.1 


0.41 


1.59 


0.24 


202 


45 


8.6 


18.6 


1.7 


0.46 


1.9 


0.29 


256 


50 


12.3 


30.8 


2.4 


0.51 


2.2 


0.34 


316 


55 


16.8 


47.6 


3.2 


0.57 


2.5 


0.39 


382 


60 


22.3 


69.7 


4.3 


0.62 


2.8 


0.44 


455 


65 


28.6 


101 


5.4 


0.67 


3.1 


0.49 


534 


70 


35.7 


150 


6.7 


0.72 


3.5 


0.54 


619 



discussed energy range. At energies higher than 15 MeV Ug detection area of LSD is 16.5 
times less than those of KII. 

Unfortunately, we are not aware of any paper in which the cross-section of the antineutrino- 
iron reaction (8) is tabulated. That is why we can not calculate non-electron antineutrino 
LSD effective area. However, cross-sections of interactions of antineutrinos with heavy 
nuclei are in general not greater than those of neutrinos (see, for example, [21]). There- 
fore we may assert that non-electron antineutrino LSD effective area was not greater than 
electron antineutrino LSD effective area: 

Sp.<S.^. (1) 

It is easily seen from the Table 6 that 

(1) KII was roughly an order of magnitude more sensitive to jy^ than to and h'rc] 

(2) electron neutrino (ve) sensitivities of LSD and KII were comparable in the energy 
range £^ ~ (25 - 45) Mev; 

(3) non-electron neutrino {vx) sensitivity of LSD was greater then those of KII for energy 
> 30 MeV. 

At this stage one may think that no contradiction between LSD and KII event numbers 
would occur if supernova neutrino flux was composed of electron and non-electron neu- 
trinos of appropriate energies. The question is weather such flux content is possible in 
principle. To answer this question one should consider neutrino emission by the collaps- 
ing supernova core and neutrino flavour transformation in the matter of the star due to 
MSW-effect. 



6 



3 Supernova neutrino emission and flavour transfor- 
mations 



There is a number of reviews on the supernova neutrino emission and flavour transforma- 
tions (see, for example, [2l], [25]). In this section we utihse the results described there in 
detail. 

Neutrinos and antineutrinos of all three flavours can be created during the collapse of 
the iron core in the center of the star. All reactions in which they are created conserve 
lepton flavour. In particular, muon neutrinos and antineutrinos are created only in pairs, 
the same is valid for tau-neutrinos and antineutrinos. Moreover, muon neutrinos are 
produced in the same reactions as tau-neutrinos. Therefore in any collapse model neutrino 
fluxes satisfy the following conditions: 

po ^ po ^ po po ^ po ^ po po ^ po (2) 

juCore 

0<^°-^^<^ (3) 

Here F^^ is a time- and energy-integrated flux of an (anti-)neutrinos, subscript specifying 
the type of (anti-)neutrinos. is a standard notation for non-electron neutrino original 
flux. Superscript "0" denotes that it is an original flux, i.e. such a flux which would reach 
the earth if there were no flavour transformations in the matter of the star. L = 52kpc is 
the distance between the supernova and the earth, and N^"^"^ is the number of electrons 
in the iron core. 

^Core ^ 26^^, (4) 

where mpe is the iron nucleus mass, and 26 is the number of electrons in the atom of iron. 

The supernova collapse begins when the iron core reaches the Chandrasekhar mass 
limit, Mch- Chandrasekhar mass varies from IAMq to 2.2Mq, Mq being the mass of the 
Sun (see, for example, review [2S] or paper [27]). Last inequality in reads 

< F° - < 0.34^^ ■ lO^W^ (5) 



Due to the MSW-effect [TT][T2] neutrinos may change their flavour while passing 
through the envelope of the star. In general neutrino and antineutrino fluxes at the 
earth Fy^^y^ are linear combinations of original fluxes [28] [29] : 



F.^=pF° +(l-p)F°, (6) 

F-y^=vFl^{\-V)Fl (7) 

Fy^ + Fy^ = (1 - v)Fl + (1 + vWl (8) 

Fp,^ + F,^ = (1 - v)Fl + (1 + V)Fl (9) 

where conversion probabilities p depend on the unknown neutrino mass hierarchy and 
neutrino mixing angle One may distinguish three extreme (with respect to 6*13) cases: 
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Table 7: Conversion probabilities p and p for different neutrino mass hierarchies and 
values of 613. 





^13 < 3 ■ 10-3 


^13 > 3 ■ 10-2 


3 • 10-3 < ^13 < 3 ■ 10-2 


Normal hierarchy 


p = sin^ 6*12 

p = COS^ 9i2 


p = 

p = COS^ 612 


0<p< sin2 012 

p = COS^ 612 


Inverted hierarchy 


p = sin^ 9i2 
p = 


p = sin^ 612 
0<p< cos2 Ou 



LN ( Large angle, Normal hierarchy), 
LI (Large angle. Inverted hierarchy), 
SA (Small angles. Any hierarchy). 

Here large angles stand for O13 > 3 ■ 10-^, while small angles stand for 6*13 < 3 • lO^^. In 
addition a case of intermediate 6*13 should be considered. Note that current experimental 
limit is ^13 < 0.17, sin^ 613 < 0.03 ^SOj. Coefficients for all this cases are given 
in Table 7 [28] [29]. 

According to PDG [S] sin^ 9^^ = 0.28. 

Two consequences follow immediately from eq.([7]) and Table 7. 

(1) Upper bound on : 

< T^F,^ < 3.5Fp^ (10) 
1 — p 

Thus the original flux was constrained by the flux at the earth Fp^, which by itself 
was severely constrained in the considered energy range by KII non-observation of elec- 
tron antineutrinos with E > 14MeV. 

(2) Upper bound on F^^ in LN and SA cases: 

F° < = 1.4Fp^, LN and SA (11) 

Thus if either LN or SA case is realised in Nature then the original flux Fp^ was 
constrained for the same reason as F°. In the LI case Fp^ is independent from F^? . 

Numerical values for both upper bounds are given in the next section (eqs. ffT^ and 

(USD). 

Values for p, p given in Table 7 are valid if the density profile of the envelope of the 
progenitor star is static during neutrino emission time interval (~ 10 seconds). This 
assumption may be violated by a shock wave, which moves from the core through the 
envelope. If it reaches the region of neutrino conversion (resonance region), it may change 
values of p and p [T3]. It is a shock wave which blows up the star envelope, therefore it 
is an inevitable feature of the supernova explosion. However, it is not clear whether the 
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first stage of the two-stage explosion should be accompanied by the shock wave powerful 
enough to reach the resonance region. 

Simple physical considerations and detailed numerical studies [IS] , [32] , [33] , [31] show 
that 

(1) shock wave reaches the resonance region and in 2-5 seconds after the onset of the 
collapse; 

(2) in the SA case the shock wave effect is negligible; 

(3) in the LN case shock wave switches p from to some non-zero value in the interval 
[0, sin^ 6*12], but does not affect p; 

(4) in the LI case shock wave switches p from to some non-zero value in the interval 
[0, cos^ 6^12], but does not affect p. 

SN1987A emerged in the Southern hemisphere, and the neutrinos had to pass through 
the earth in order to reach LSD and KII. The earth matter effect for supernova neutrinos 
is described in detail for example in [33|. Not going into details we note that in general 
it could change p and p by 5%- 30%, but in the LN case it could not affect p, and in the 
LI case - p. As a result the earth matter effect proved not to influence significantly the 
results of the statistical analysis described in the next section. 

4 Statistical analysis and results 

Energy distribution of supernova neutrinos is thought to be some modification of Fermi- 
Dirac distribution (see, for example, [22])- Unfortunately, it is absolutely impossible to 
reconstruct the shape of the distribution from the LSD neutrino signal because as it 
was stressed above the measured event energy carried almost no information about the 
incident neutrino energy, to say nothing about small number of events. Therefore in our 
statistical analysis we consider a set of " monoenergetic" hypothesises. To be more precise, 
we numerically investigate the probability P{E, Fj^^, F^^, F^^ ) that monoenergetic (with 
energy E) fluxes F^? and F° after the flavour transformation according to eqs.(l6])-(l9]) 
produced 

(1) not less than 5 events in LSD, 

(2) not greater than 2 events in KII, no one from which being a ue elastic scattering with 
Ee > 14 MeV. 

Poisson statistics is applied. Energy runs from 15 to 70 MeV with the step 5 MeV. Three 
combinations of neutrino mass hierarchy and mixing angle ^13, described in the previous 
section, are implemented. 

First we maximize P{E, F^^, F^^, F^) with respect to F°^, fixing each of the other two 
original fluxes to be one of the values from the set {(0.5, 1, 1.5, 2) x lO^cm"^}. We find 
that 

1. the maximal value of probability P does not exceed 1% for values of F° greater than 
2 ■ lO^cm"^ for all combinations of neutrino parameters (including extreme cases LN , 
LI , SA ): 

F° < 2 ■ lO^cm-^, 99%CL; (12) 

2. the maximal value of probability P is less than 1% for values of F^^ greater than 
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Table 8: Maximal probability, lower and upper bounds and central value for the effective 
original flux in the LI and SA cases. Fluxes are measured in 10^'' cm^. See text for further 
explanations. 



E^{MeV) 


p 

^ max 


lower 


central 


po 

upper 


30 


0.030 


2.2 


3.4 


5.0 


35 


0.037 


1.2 


2.1 


3.3 


40 


0.032 


0.77 


1.2 


1.8 


45 


0.024 


0.57 


0.75 


0.97 



lO^cm ^ for the LN and SA cases: 

< 10^cm-^ 99%CL; (13) 
Combining this with inequality ([5]) one obtains an upper bound on : 

po < 0.34^^ ■ 10l°cm-^ LN or SA . (14) 

Simple estimates of the mean number of events in the detectors, as well as final re- 
sults, confirm that the value 10®cm~^ is negligible compared to the original neutrino flux 
sufficient to explain LSD signal, which is of order of 10^°cm~^. Thus, in the subsequent 
analysis we may safely put = and, in the LN and SA cases, F^^ = 0. In order to 
consider the LI case on equal footing with the LN and SA cases, we note that according 
to (HD 

P{E, Fl F° , 0) < PiE, Fl + F° , 0, 0). (15) 

This allows to consider P{E, + Fp^, 0, 0) as an upper bound on P{E, F°^, F^? , 0) in the 
LI case. 

After this simpliflcations we maximize probability F(F, F°, 0, 0) with respect to effec- 
tive original flux F°, which stands for F^ in the LN and SA cases and for F° +F^, in the 
LI case. When the maximal probability Fmax is greater than 2% we find lower and upper 
bounds for F°, F°ppgj, and Fj^^^j. correspondingly, solving equation F(F, F°, 0, 0) = 0.02. 
The result for the LI and SA cases is presented in Table 8, and for the LN case - in 
Table 9. Only those energies, for which Fmax > 2%, are included in the tables. 

The following conclusions may be derived from this tables 

SA case. Those values of Fp^, which agree with the data (see Table 8), do not fill 
inequality ffT^ . The shock wave effect, as noted in the previous section, is negligible. 

LI case. Electron neutrinos and antineutrinos in the energy range (30-45) MeV with 
the total flux of order of 10^*^ cm^ flt the data. Inequality ( |T^ does not hold for this 
case, the zero value of p = being the reason for this. However, a powerful shock wave 
could switch p to non-zero value, thus making inequality (1141) obligatory for this case also. 
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Table 9: Maximal probability, lower and upper bounds and central value for the effective 
original flux in the LN case. Fluxes are measured in 10^° cm^. The lower block of the 
table correspond to the flux allowed by inequality IHM with Mcore = 2Mq. See text for 
further explanations. 



E^{MeV) 


p 

max 


lower 


central 


Upper 


25 


0.029 


9.6 


14.7 


21.2 


30 


0.075 


4.1 


9.7 


19.5 


35 


0.13 


2.3 


6.7 


16.3 


40 


0.19 


1.4 


4.8 


13.5 


45 


0.25 


0.91 


3.5 


10.9 


50 


0.31 


0.64 


2.7 


8.7 


55 


0.36 


0.46 


2.1 


6.9 


60 


0.40 


0.35 


1.64 


5.6 


65 


0.44 


0.27 


1.34 


4.6 


70 


0.48 


0.22 


1.11 


3.8 



Therefore shock wave could spoil the agreement of the theory with the data. 

LN case. Taking into account inequality ffT^ with Mcore = SM©, one may conclude 
that electron neutrinos with energies greater than 60-65 MeV with the total flux of order 
of 10^° cm^ fit the data. A powerful shock wave could switch p from to some nonzero 
value, thus effectively switching the LN case to some modification of the SA case. This 
could spoil the agreement of the theory with the data. 

Maximal value of probability in Table 8 does not exceed 4%. Table 9 and inequality 
(fT4l) imply unusually high neutrino energies. Evidently, agreement of the theory with the 
data is far from being perfect. However, various possible sources of inaccuracy in the 
analysis may be suspected, for example the following. 

(1) Cross section of neutrino- nucleus interactions used above are obtained mainly as a 
result of calculations, not measurements, and may appear to be inaccurate. Thus LSD 
effective area may appear to be greater, and KIT effective area - smaller than in Table 6. 

(2) The value of the efficiency of LSD neutrino event reconstruction which we utilised 
according to [lOj, 53%, is claimed to be an estimate. If it is greater than 53%, then the 
LSD effective area increases. 

To show that our results are stable under the reasonable variation of the input, we 
repeated calculations with the increased LSD effective area (by factor 1.3) and decreased 
KII effective area (by factor 1.4). We believe that such a procedure accounts for possible 
input variations in an optimistic (with respect to the concordance of the LSD signal with 
the double-burst hypothesis) way. The results are presented in Tables 10 and 11. In this 
tables two pairs of upper and lower bounds are presented, one of which corresponding 
to P = 2%, and another - to P = 5%. It is clear from this tables that most of our 
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Table 10: Maximal probability, lower and upper bounds and central value for the effective 
original flux in the LI and SA cases, in case of 1.3 times increased LSD and 1.4 times 
decreased KII effective areas. Fluxes are measured in lO^'' cm^. The lower block of the 
table correspond to the flux allowed by inequality ( |T4l) with Mcore = 2Mq. See text for 
further explanations. 



E^{MeV) 


pU 

max 


lower 2% 


T— ifl 

lower 5% 


central 


upper 5% 


T-rfl 

-'^ upper 2% 


25 


0.031 


3.4 




5.3 




7.8 


30 


0.066 


1.47 


2.3 


3.3 


4.5 


6.2 


35 


0.089 


0.82 


1.20 


2.0 


3.2 


4.2 


40 


0.092 


0.51 


0.74 


1.27 


2.0 


2.6 


45 


0.082 


0.34 


0.50 


0.81 


1.2 


1.6 


50 


0.070 


0.24 


0.37 


0.54 


0.75 


1.01 


55 


0.058 


0.18 


0.29 


0.37 


0.47 


0.67 


60 


0.049 


0.14 




0.27 




0.46 


65 


0.039 


0.11 




0.20 




0.31 


70 


0.026 


0.10 




0.14 




0.19 



Table 11: Maximal probability, lower and upper bounds and central value for the effective 
original flux in the LN case, in case of 1.3 times increased LSD and 1.4 times decreased KII 
effective areas. Fluxes are measured in 10^° cm^. The lower block of the table correspond 
to the flux allowed by inequality ( fT^ with Mcore = 2Mq. See text for further explanations. 



E^{MeV) 


p 

max 


^ lower 2% 


^ lower 5% 


central 


upper 5% 


-'^ upper 2% 


25 


0.055 


6.4 


10.9 


13.2 


15.9 


23.8 


30 


0.12 


3.0 


4.2 


8.5 


15.5 


19.8 


35 


0.19 


1.7 


2.3 


5.8 


12.9 


16.0 


40 


0.26 


1.0 


1.4 


4.0 


10.5 


12.6 


45 


0.33 


0.69 


0.92 


2.9 


8.3 


9.6 


50 


0.38 


0.48 


0.64 


2.2 


6.5 


7.4 


55 


0.44 


0.35 


0.47 


1.70 


5.1 


5.8 


60 


0.48 


0.27 


0.35 


1.35 


4.1 


4.6 


65 


0.52 


0.21 


0.28 


1.10 


3.3 


3.7 


70 


0.55 


0.17 


0.22 


0.91 


2.8 


3.1 



12 



conclusions remain unchanged. Only the status of the SA case changes from "excluded" 
to "disfavoured". Note that in all cases but one, the LI case without shock wave, large 
mass of the iron core M'-"°'''^ is favoured. 

5 Conclusions 

The following conclusions concerning the first stage of the two-stage SN1987A explosion 
models are obtained. 

(1) In the case of small mixing angle ^13, ^13 < 0.003, such models are disfavoured by the 
data, independently of the neutrino mass hierarchy. 

(2) In any model non-electron neutrino and antineutrino production had to be severely 
suppressed during the first stage of the collapse, independently of the neutrino mass hi- 
erarchy and mixing angle ^13 : 

F° < 10«cm-2. (16) 

This means that at the first stage of the collapse there was no thermal equilibrium, even 
rough. 

(3) In the case of normal mass hierarchy and large mixing angle ^13, ^13 > 0.03, in 
order to explain the data one should imply emission of very energetic {E > 60 MeV) 
electron neutrinos, z/g, at the first stage of the explosion; at the same time the suppression 
of Ue production should be assumed: 

F° ^ (0.3 - 0.5) • 10^°cm-2, F^^ < lO^^cm"^. (17) 

In addition, large values of the collapsing core mass, Mcore ^ 2M0, are necessary. A 
powerful shock wave could further complicate the agreement of the data with the theory. 

(4) In the case of inverted mass hierarchy and large mixing angle ^13, ^13 > 0.03, the 
data can be explained by the moderate energy (30 MeV< £^ < 45 MeV) electron neu- 
trino and antineutrino emission at the first stage of the explosion with fiuxes of order of 
lO^'^cm"^. A powerful shock wave could worsen the agreement of the data with the theory. 

The conclusions of the current work incorporate those obtained in the previous in- 
vestigation [lAl of the author. In particular, following ^lOj and ^/Tj, in [H] the energy 
of the neutrinos responsible for the LSD signal was supposed to be in the 30-50 MeV 
range; under this assumption it was concluded that the case of normal mass hierarchy 
and large mixing angle ^13 is not compatible with the data. This accords with our current 
conclusion (3). 
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